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We consider the Kondo lattice model in two dimensions 
at half filling. In addition to the fermionic hopping integral t 
and the superexchange coupling J the role of a Coulomb re- 
pulsion U in the conduction band is investigated. We find the 
model to display a magnetic order-disorder transition in the 
U-J plane with a critical value of Jc which is decreasing as 
a function of U. The single particle spectral function A{k,u}) 
is computed across this transition. For all values of J > 0, 
and apart from shadow features present in the ordered state, 
A{k,u!) remains insensitive to the magnetic phase transition 
with the first low-energy hole states residing at momenta 
k = (±7r,±7r). As J ^ the model maps onto the Hub- 
bard Hamiltonian. Only in this limit, the low-energy spectral 
weight at fc = (±7r,±7r) vanishes with first electron removal- 
states emerging at wave vectors on the magnetic Brillouin 
zone boundary. Thus, we conclude that (i) the local screen- 
ing of impurity spins determines the low energy behavior of 
the spectral function and (ii) one cannot deform continuously 
the spectral function of the Mott-Hubbard insulator at J = 
to that of the Kondo insulator at J > Jc. Our results are 
based on both, T = Quantum Monte-Carlo simulations and 
a bond-operator mean-field theory. 
PACS numbers: 71.27.-Ha, 71.10.-w, 71.10.Fd 



I. INTRODUCTION 

Starting from the seminal work of Brinknian and 
Rice [|l[ the single-hole dynamics in correlated insulators 
has remained to be an intriguing issue with many open 
questions yet to be clarified. In this respect Kondo lat- 
tice systems at half filling provide for a case in which 
the single particle dynamics can be studied continuously 
across genuinely distinct correlation induced insulating 
phases, i.e., Mott-Hubbard and magnetic insulators as 
well as Kondo insulators. Of particular interest in this 
situation is to understand i) which properties of the cor- 
related insulator, i.e., long range magnetic order or local 
effects determine the functional form of the quasiparticle 
dispersion relation and more specifically ii) if it is pos- 
sible to continuously deform the spectral function of the 
Mott Hubbard insulator to that of the Kondo insulator. 
In order to answer these questions, it is the purpose of 
this paper to consider a Kondo lattice model (UKLM) on 
a two-dimensional square lattice with an additional local 
Coulomb repulsion U between the conduction electrons 
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The unit cell, i, contains a localized orbital and an 
extended conduction band state. In the Kondo limit, 
charge fluctuations on the localized orbital are sup- 
pressed with only the spin degrees of freedom remaining. 
Si = ^s s' fi ^s,s'fi //2, where a are Pauli spin-1/2 

matrices and /J are fermionic operators which satisfy 
the constraint y^^flf- = 1. Conduction band elec- 
trons of spin z-component a are created by cl where 
n^ ~ cl c- is the conduction band density for spin 
z-component a . The extended orbitals overlap to form a 
band with a dispersion e(A:) = — 2t(cos(fc2;) -f cos(fcj,)) as- 
suming a nearest- neighbor (NN) hopping integral i. The 
Coulomb repulsion V is taken into account by the Hub- 
bard interaction term. 

At half-filling and for the particular conduction band 
structure chosen the UKLM is an insulator for all values 
of U and J [|-|. Specifically, as J ^ the UKLM 
maps onto the Hubbard model. The latter is in a Mott 
insulating phase, which however due to nesting on the 
square lattice is masked by an antiferromagnetic spin 
density wave with no spin gap and a charge gap oc [/ 
in the strong coupling limit. As J jlJ —* oo the Hub- 
bard repulsion can be neglected relative to the exchange 
scattering and the model maps onto the pure Kondo lat- 
tice model (KLM) with J jt > 1. In this limit, the 
ground state is a Kondo insulator with spin (A^) and 
charge (Ac) gaps satisfying Ac > A^ Q. In both of the 
aforementioned limiting cases, the single particle spec- 
tral function displays very different behavior. The lim- 
iting ground state for J jt -^ c», i.e. the Kondo insu- 
lator, is a direct product of Kondo singlets on the /-c 
bonds of the unit cell. Adding a hole into the conduction 
band will break a singlet and leads to a hole dispersion 
relation e(/c) — iJ/A -f t(cos(/ca;) -f- cos(A;j,)) in first or- 
der perturbation theory in t/ J |^]. Hence the first elec- 
tron removal-states ('Fermi points') occur at wave vec- 
tors k = (±7r, ±7r) within the Brillouin zone. At present 
the precise form the the single particle spectral function 
for the Mott insulating state is still unknown. Yet, var- 
ious numerical and analytical approaches confirm that 



the first electron removal-states are located at fc-points 
on the boundary of the magnetic Brillouin zone, i.e., at 
k = (0,±7r),(±7r,0) and k = (±7r/2, ±7r/2). In order 
to shed light onto this situation we can fix U and, as 
a function of J/t, drive the system through a magnetic 
quantum phase transition at J = Jc{U) from the Kondo 
insulator for J/t >> 1 into the Mott insulating state for 
J ^ 0. Along this path we compute the spectral function 
^(^,0;), both, exactly by using Quantum Monte Carlo 
methods and approximately using a bond-operator mean 
field theory. Based on our findings we will argue that the 
low energy features of the spectral function are insensitive 
to the quantum phase transition. Speaking differently, 
the low energy hole-states are found at fc = (±7r, ±7r) for 
all values of J > 0. It is only at J = that the spectral 
weight of the low energy feature at fc = (±7r, ±7r) van- 
ishes to produce the single-hole dispersion relation of the 
Hubbard model. Thus our main results are (i) that the 
local screening of the /-spins dominates the low energy 
features of the spectral function and (ii) that there is no 
continuous path from the Kondo to the Mott-Hubbard 
insulator in this specific model. 

The paper is organized as follows. In the next section, 
we briefly outline the Quantum Monte-Carlo (QMC) 
method as well as the bond-operator-mean field theory. 
In section [II, we first discuss the magnetic phase dia- 
gram of the UKLM model. We map out the critical line 
in the U-J plane to show that Jc is a monotonically de- 
creasing function of U. Comparison of the QMC and 
mean-field results prove to be very satisfactory. Hav- 
ing determined Jc as a function of U/t we then focus on 
the single particle spectral function and detail the afore- 
mentioned results (i) and (ii). Section IV is devoted to 
conclusions. 



the reader to Refs. [Q-g] for a detailed description of the 
method. 



B. Mean field theory 

For an approximate description of the UKLM we apply 
a mean field theory similar to the one proposed for the 
pure KLM in m, where further details can be found. 
We represent the local Hilbert space consisting of one / 
electron and additionally up to two itinerant electrons by 
applying the following operators onto the vacuum |0) of 
an empty site 

^^|o> = ^(4/! + /tM)Io) 

4|0) = ^(c|/|-ct/t)|0) 
tt|0) = -j=(ct/t+ct/t)|0) 



4|0) = ^(c|/|+ct/t)|0) 
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4lo> = /i|o) 

bUo) = c\cin\o) 



(2) 



The s and t operators are equivalent to the so-called bond 
operators of ||8|| and are assumed to obey bosonic commu- 
tation relations. The fermionic operators a and b have 
been introduced first in |^,|0[ and label states with one or 
three electrons per site. In order to suppress unphysical 
states a constraint of no double occupancy 
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II. METHODS 

A. Quantum Monte Carlo 

We have used the projector auxiliary field Quantum- 
Monte-Carlo (PQMC) method to investigate the UKLM 
model. This method is based on the projection equation 
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where |^t) is required to be non-orthogonal to the 
ground state |^o)- Details of how to implement the 
PQMC method without generating a sign problem in 
the particle-hole symmetric case, i.e. at half-filling, have 
been introduced and extensively described for the KLM 
in Refs. HJ^]. We have also used a new and efficient 
method pI to calculate imaginary time displaced Greens 
functions within the PQMC formalism. From the tech- 
nical point of view no complications arise when intro- 
ducing the Hubbard term into the KLM. Hence we refer 



has to be fulfilled. 

Rewriting the UKLM in terms of (Q) leads to a strongly 
correlated boson-fermion model which cannot be solved 
exactly. To proceed we use a mean-field approach which 
incorporates both, an antiferromagnetic state and a spin- 
singlet regime. The latter regime can be expressed by a 
condensate of the singlets [g| 

(^,> = (4) = '■ (4) 

In the antiferromagnetic phase the dimers condense into 
a linear combination of the singlet and one of the triplets 
M^ implying that 
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Here ( — 1)^ is a shorthand for '+!(— 1)' on the magnetic 
A(B) sublattice. 

Inserting ^ into the UKLM and using the mean field 
approximation (y,|5|) we obtain 
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where Pf(x) = 1(~1) and we have introduced a chemi- 
cal potential A to set the global particle density and a 
local Lagrange multiplier /^^ in order to to enforce the 
constraint (0). In the remainder of this work we assume 
/ii to be site independent, i.e. /i; — fi. The difference 
between (R) and the mean-field Hamiltonian for the pure 
KLM |7| resides in the last line of (|^) which accounts 
for a suppression of doubly occupied conduction electron 
orbitals. 

Diagonalizing (R) leads to 4 bands LJi2{k) = X±Ei{k) 

and uj3,4{k) = XztE2{k) which are twofold degenerate by 
spin-z quantum number 
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Here e^: 



'2/^X](i=i cos fed . Note, that the dispersions in 
(^ do depend on U, as the order parameters s, ?7i and /i 
are functions of t/. At half filling the lower(upper) two 
bands, i.e 102,4 (fj^i.s), are completely filled(empty). 

Evaluating the ground state energy and using the sta- 
tionarity conditions dE/ds = 0, dE/dm = 0, and 
dE/dfi = the mean-field equations in the magnetic 
phase {m ^ 0) read 
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FIG. 1. Mean-field order parameter s, m and the staggered 
magnetizations Mc and Mf as function of (a) J/t for U/t — 4 
and (b) U/t for J = t. 
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Fig. |l| shows numerical solutions of (g) and ( [lO|) as 
a function of J and [/. In both cases we find a second 
order phase transition between the antiferromagnetically 
ordered and the Kondo phase. Fig. 1 also shows the 
staggered magnetizations Mc(f) of the c(/) electron 
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Using (0) we may express the spectral function 

Ac{k, uj) of the conduction electron cg via a multi-particle 
correlation function of the s, t, a and h operators. On 
the mean- field level however, this simplifies into a lin- 
ear combination of one-particle propagators of the a and 
b fermions only, involving both, diagonal as well as off- 
diagonal contributions. After some algebra we get 

Ac{k,z)^ (12) 
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III. RESULTS 



where Et = 2(E l± E i). For the disordered Kondo- 
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singlet phase (m = 0) we get 
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FIG. 2. Phase diagram of the Kondo-Hubbard modeL 
Sohd hne: QMC, dot-dashed hne: MF, dashed hne: Spin 
Hamiltonian Hspin Eq. (|l3|) . Circles show parameter values, 
where the spectral function has been evaluated. 



A. Magnetic phase diagram 

We start the discussion of our results with the mag- 
netic phase diagram. At f/ = the UKLM maps 
onto the KLM. In the latter, the competition between 
the Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion and the Kondo screening leads to a quantum phase 
transition between an ordered magnetic state and the 
disordered singlet phase at Jc/t ^ 1.5 jl^,^,^. For 
U/t -^ oo double occupancy of the conduction electron 
sites is suppressed and the model maps onto a pure spin 
Hamiltonian of the form: 
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with J|| = At^/U. Hence in the limit U — > oo, Jc vanishes 
and the ground state is a product of singlets on the f-c 
bonds. 

We have determined Jc as a function of the Hubbard 
repulsion U both, on the mean-field (MF) level and with 
the QMC method. Within MF theory the staggered mag- 
netization is given by Eq. (fl), while from the QMC 
method it is determined using the static spin-spin corre- 
lation function 
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a — c{f) labels conduction(f ) electron spins S^ and the 

total spin SiP* is given by 5"^°* = SS + sL The staggered 
moment is extracted from finite size extrapolation 
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FIG. 3. Staggered magnetization for a fixed J = O.i 



Fig. g depicts the phase diagram as a function of J/t 
for finite U/t. The solid line refers to QMC results, the 
dashed-dotted line shows the mean- field results. As an- 
ticipated already by the preceding discussion of the lim- 
iting points C/ — > cx) and U = 0, the critical value J is 
a monotonically decreasing function of U. This can be 
understood as the Hubbard interaction tends to localize 
the conduction electrons leading to an effective reduction 
of the hopping amplitude. Hence, the formation of local 
singlets is favored. The above spin Hamiltonian ([l3|) has 
been analyzed by Matsuhita and collaborators |13[ who 
find a phase transition between a spin liquid and anti- 
ferromagnetically ordered phase at (J||/J)c = 0.71. This 

leads to Uc — T^ry. the dashed hne in fig. |[ 



m consis- 



0.71J' 

fence with the results for the Kondo-Hubbard model in 
the large U/t limit. 

Finally Fig. plots the staggered magnetization from 
a QMC scan at fixed J. The broken-symmetry ground 
state satisfies m*°* = m^ — m^. In QMC 77i*°* was cal- 
culated independently and up to [/ < 4 the above re- 
lation is fulfilled within the errorbars. With increasing 
U conduction electrons get more and more localized and 
their local moment grows until it reaches the maximum of 

/ (S*?)^ ) — 3/4 in the strong coupling region U > 8. The 

staggered moments in the small f/ < 4 region are well 
understood within a Neel picture of almost fully ordered 
/-spins where the small local moment of a conduction 
electron is anti-parallel to the impurity spin. For larger 
values of U dimerization becomes important which sup- 
presses both m'^-^ . 



B. Single particle spectral function. 

To study the single-hole dynamics we analyze the spec- 
tral function A{k,uj), both using the MF expression (|l^ 
as well as results from the QMC. Within the QMC ap- 
proach we first evaluate the imaginary time Greens func- 
tion 



where Q = (tt, tt) and iV is the number of unit cells. 
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FIG. 4. Single particle spectral function for pure Hubbard 
model at [/ = 4, (7 = 6 and (7 = 8. For the QMC data (solid 
lines) we have normalized the maximum peak heights to unity. 
The numbers on the left-hand side of the figures correspond 
to the normalization factor. The vertical bars show the MF 
dispersion relation. 



from which A{k,uj) is extracted using the maximum en- 
tropy (ME) method Q. 

We begin with the pure Hubbard model. In Fig. ^ 
we plot Aik^bj) as obtained from QMC as well as the 
MF dispersion as a function of U jt. While the compari- 
son of the QMC with the MF dispersion is favorable one 
has to realize that the MF approach overestimates the 
quasiparticle gap. Therefore the MF band structure in 
these figures results from taking only s and m as obtained 
from the self-consistency equations (g) however adjusting 
^ such as to obtain the QMC gap at fc = (f , f ). At weak 
coupling U /t <C 1 we find that the overall form of the 
low-energy dispersion is well reproduced by a functional 
form ±(A^ -I- eilz)^/"^ which is consistent with that in a 
spin density wave (SDW) state. Exactly this dispersion 
emerges also from the bond-operator MF theory at J = 
where E^ ^ reduces to the SDW dispersion and the spec- 
tral weight of excitations with the dispersion £„ t van- 
ishes. For J = the condensate densities for the triplet 
and singlet are identical, i.e. m = s, which is equivalent 
to a Neel state of the /-spins. 

At strong coupling U /t ^ 1 the Hubbard model maps 
approximately onto the t-J\\ model with J|| = AP/U. 
Monte Carlo results for the latter model at J||/i < 1 
show the existence of a quasiparticle band of width ^ Jn 
||l5| . This should be compared to an identical spectral 
feature which can be observed in our QMC data for the 
Hubbard model upon enhancing U/t in fig. yb,c) (see 
also [0). Especially along the line from k — (0,7r) to 
k = (0; 0), this narrow quasi particle band is clearly visi- 
ble. In principle one should observe a similar band along 
(0, 7r)-(7r, tt), however, due to small spectral weight in 
this region, we are unable to resolve this feature. Of par- 
ticular importance is, that for the parameters we have 
investigated, the momenta of the dominant lowest en- 
ergy hole-states for the Hubbard model are found on the 
boundary of the magnetic Brillouin zone. For the cal- 
culations presented in this work we have been unable to 
resolve an energy difference between the (7r/2,7r/2) and 
(0,7r) points. 

Next we turn to the UKLM at finite J. In fig. | 
we show a scan of QMC spectral functions and the MF 
dispersion ranging from the Kondo phase for J/t =1.5 
and U/t = 4 to the antiferromagnetic phase at J/t = 
0.4, 0.6 and U/t = 4. As for the pure Hubbard model the 
QMC and MF results are reasonably consistent. From 
the perturbative argument for J/t >> 1, given in section 
ffl, we expect the momenta of the dominant lowest energy 
hole-states to occur at fc = (±7r, ±7r). As can be seen 
from fig. I^a) this is in consistent, both with the QMC 
as well as with the MF results. Moreover the QMC and 
MF dispersion agree very well. 

Lowering J as in fig. |^a)-c) reveals the evolution of the 
spectral density on going from the Kondo to the antifer- 
romagnetically ordered phase. In fact, as J approaches 
zero additional bands with a dispersion similar to the 
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FIG. 5. Single particle spectra for the Kondo phase (a) and 
in the antiferromagnetic phase (b,c). Vertical bars show the 
MF dispersion. 



pure Hubbard case, i.e. ^), develop. Yet, in the anti- 
ferromagnetically ordered phase, but for a finite J the 
lowest energy hole-states are still Kondo-like, i.e. they 
occur at fc = (±7r, ±7r) as can be seen in fig. 0b), c). How- 
ever, the weight of this excitation decreases continuously 
with decreasing J and vanishes at J = 0. The weight 
of the Hubbard-like band at fc = (tt, tt) increases from 
zero in the spin singlet phase to its maximum value at 
J = 0. Therefore we can interpret the change in the spec- 
tral function with decreasing J as a continuous transfer 
of weight from Kondo-like to Hubbard-like bands. This 
shift of spectral weight renders the J = point singular 
since there is a sudden change of the wave vector which 
dominates the low energy hole dynamics. 




FIG. 6. Relative spectral weight of Hubbard-band vs. 
Kondo-band. Thick lines: MF. Thin lines: QMC. 

These findings are corroborated by our MF results. In 
fig. ||the relative weight rH(K) — Zh{k)/{Zh + Zk) of 
the lower(upper) Hubbard (Kondo)-like band at momen- 
tum k = (tt, tt) as obtained from integrating Ac{k,Lo) is 
depicted. In the QMC approach Zk results from fitting 
the long-time tail of the Greens function at fc = (tt, tt) 
to the form Zp. "" where A^p corresponds to the quasi- 
particle gap. In turn Zh is obtained from the sum rule 
Zk + Zh = Tm{k) assuming a two-pole structure. Both, 
the QMC and the MF approximation display the same 
overall trend: at J/t = the total weight is in the 
Hubbard-like band while with increasing J it becomes 
distributed into both bands. In the Kondo phase the 
Hubbard band disappears completely. In addition fig. o 
shows, that the MF approximation underestimates the 
spectral weight in the Hubbard-like band. 

To compare the momentum dependence of the spectral 
weight as obtained from the MF theory with that of the 
QMC fig. depicts Ac{k,uj) from (|l|) for U/t = 4 and 
J/t = 0.4. For visualization purpose, we have smeared 
the delta-function like MF-spectrum by a finite imagi- 
nary part S = 0.03. The weight of these delta-peaks 
strongly varies as a function of k having its maximum 
around (0, 0) and a very small value in the vicinity of 
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FIG. 7. Mean-field spectral lunction Eq. [13 in the antifer- 
romagnetic phase. 



(tt, tt). Again, this is consistent with the QMC data of 
fig. He). Obviously, since the imaginary part of the self 
energy vanishes in the MF approximation, the broaden- 
ing of the QMC spectral function is not reproduced. Note 
however, that on the QMC side pinning down the details 
of the line shape is extremely challenging. 



erate ground state. In turn, the evolution of the spectral 
function is discontinuous in as such that at J = there 
is a sudden change of the wave vector which dominates 
the low energy hole dynamics. In this sense the model 
shows no continuous path from the Kondo insulator to 
the Mott insulator. 

The singularity of the UKLM at J = may be allevi- 
ated by including an antiferromagnetic coupling between 
the localized /-spins. In the large U/t limit this leads 
to a bilayer spin model which has been considered by 
Vojta and Becker p^. The authors arrive at a similar 
conclusion namely that hole dynamics are governed by lo- 
cal spin environment. Numerical work on this modified 
model is presently under progress. 

Given our results it is very tempting to speculate on 
the effects of doping with a finite density of holes rih away 
from half filling. In the limit J/i — » oo the Kondo lat- 
tice model can be mapped onto a Hubbard model with 
a Coulomb repulsion m ^ oo and a particle density Uh 
p8[ . In this low-density limit single particle renormaliza- 
tions iQ may be neglected which suggests that doping 
the UKLM can be understood approximately within a 
rigid-band picture. From this we would conclude that off 
half filling the UKLM displays a Fermi surface centered 
around k = (±7r, ±7r) for all values of U and J > 0. 



IV. CONCLUSION 



We have considered the single-hole dynamics in the 
Kondo-Hubbard model using both, QMC methods and 
a bond-operator mean-field approximation. Both ap- 
proaches allow for similar conclusions. The UKLM shows 
a magnetic order-disorder transition. At C/ = this 
transition is triggered by the competition between the 
RKKY interaction and the Kondo screening and occurs 
at Jc/t ~ 1.5. In the large U/t limit the model maps 
onto a bilayer spin-model and Jc scales to zero. Our 
results show that both limiting cases are linked continu- 
ously and that Jj, is a monotonically decreasing function 
of U . Hence as far as Jc is concerned the dominant effect 
of the Hubbard interaction U is to localize the conduction 
electrons which favors screening of the localized spins. 

The single particle spectral function was shown to be 
insensitive to the magnetic phase transition. Irrespective 
of U and J > the dominant low energy hole-states are 
found at the momenta k = (±7r,±7r). These excitations 
originate from the screening of the magnetic impurities 
and hence are local. In the ordered phase pronounced 
shadow features can be observed. As J ^- 0, the spec- 
tral weight in the Kondo-like low energy band in the 
vicinity of fc = (iTr, ±7r) vanishes and is transfered to 
higher energy Hubbard-like bands. In the {U, J) -plane 
the Hubbard-line, i.e. J = 0, is singular since the local- 
ized spins decouple and lead to a macroscopically degen- 
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